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Abstract 


Ship  roll  motions  in  waves  can  be  significant  due  to  small  roll  damping  and  the 
proximity  of  ship  natural  roll  frequency  to  encountered  wave  frequencies.  Roll  motions 
can  be  reduced  by  various  methods,  including  bilge  keels  and  active  roll  stabilization 
using  rudders  or  dedicated  stabilizer  fins.  U-tube  tanks  and  flume  tanks  can  also 
be  used  to  reduce  roll  motions,  and  have  the  advantages  of  being  effective  at  zero 
ship  speed  and  being  protected  from  environmental  hazards  such  as  ice.  This  report 
describes  the  implementation  of  a  U-tube  tank  model  for  the  ShipMoSD  ship  motion 
library.  The  U-tube  tank  model  has  been  implemented  for  computations  in  both  the 
frequency  and  time  domains.  Example  computations  for  a  generic  frigate  demonstrate 
the  reduction  of  roll  motions  using  a  U-tube  tank.  The  relative  effectiveness  of  a  U- 
tube  tank  will  be  greatest  when  a  ship  has  small  roll  damping,  such  as  a  vessel  with 
small  or  no  bilge  keels. 


Resume 


Le  mouvement  de  roulis  des  navires  dans  les  vagues  pent  etre  important  a  cause  du 
faible  amortissement  du  roulis  et  de  la  proximite  de  la  frequence  naturelle  de  roulis 
du  navire  par  rapport  a  la  frequence  de  roulis  des  vagues  rencontrees.  Le  mouvement 
de  roulis  peut  etre  reduit  a  I’aide  de  divers  moyens,  y  compris  I’installation  de  quilles 
de  roulis  et  de  dispositifs  de  stabilisation  actifs,  comme  des  gouvernails  ou  des  derives 
de  stabilisation.  Les  citernes  a  tube  en  U  et  les  citernes  antiroulis  peuvent  egalement 
etre  utilisees  pour  reduire  le  mouvement  de  roulis  ;  elles  ont  I’avantage  d’etre  efficaces 
meme  lorsque  le  navire  est  a  I’arret  et  sont  protegees  contre  les  dangers  associes  a 
certaines  conditions  environnementales,  comme  la  glace.  Le  present  rapport  decrit  la 
mise  en  oeuvre  d’un  modele  de  citerne  a  tube  en  U  destine  a  enrichir  la  bibliotheque  de 
donnees  sur  le  mouvement  des  navires  du  logiciel  ShipMoSD.  Le  modele  de  citerne  a 
tube  en  U  a  ete  mis  en  oeuvre  pour  effectuer  des  calculs  dans  le  domaine  temporel  et  le 
domaine  frequentiel.  Des  exemples  de  calculs  pour  une  fregate  generique  demontrent 
la  reduction  du  roulis  lorsque  I’on  utilise  une  citerne  a  tube  en  U.  L’efficacite  relative 
d’une  citerne  a  tube  en  U  sera  maximale  lorsque  le  navire  amortit  faiblement  le  roulis, 
comme  dans  le  cas  des  navires  possedant  une  quille  de  roulis  de  petite  taille  ou  n’en 
possedant  aucune. 
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Executive  summary 


Modelling  of  U-tube  Tanks  for  ShipMoSD  Ship  Motion 
Predictions 

Kevin  McTaggart;  DRDC  Atlantic  ECR  2011-300;  Defence  Research  and 
Development  Canada  -  Atlantic;  January  201 2. 

Introduction:  Ship  roll  motions  in  waves  can  be  significant  due  to  small  roll  damping 
and  the  proximity  of  ship  natural  roll  frequency  to  encountered  wave  frequencies. 
Roll  motions  can  be  reduced  by  various  methods,  including  bilge  keels  and  active  roll 
stabilization  using  rudders  or  dedicated  stabilizer  fins.  U-tube  tanks  and  flume  tanks 
can  also  be  used  to  reduce  roll  motions,  and  have  the  advantages  of  being  effective 
at  zero  ship  speed  and  being  protected  from  environmental  hazards  such  as  ice. 

Principal  Results:  This  report  describes  the  implementation  of  a  U-tube  tank 
model  for  the  ShipMoSD  ship  motion  library.  The  U-tube  tank  model  has  been 
implemented  for  computations  in  both  the  frequency  and  time  domains.  Example 
computations  for  a  generic  frigate  demonstrate  the  reduction  of  roll  motions  using  a 
U-tube  tank.  The  U-tube  tank  capability  has  been  integrated  into  ShipMoSD  Version 
S.O. 

Significance  of  Results:  The  ShipMoSD  ship  motion  library  can  now  be  used  to 
assess  the  influence  of  U-tube  tanks  on  ship  motions.  Parametric  studies  can  be 
performed,  including  investigation  of  the  influence  of  tank  dimensions  and  associated 
natural  frequency  on  ship  roll  motions.  Frequency  domain  predictions  are  computa¬ 
tionally  efficient  and  provide  data  that  is  readily  usable.  Time  domain  predictions  are 
useful  for  examining  nonlinear  phenomena,  such  as  when  U-tube  tank  fluid  motions 
and  effectiveness  are  limited  by  the  size  of  the  tank. 

Future  Plans:  The  ShipMoSD  ship  motion  library  will  be  actively  maintained  and 
supported. 
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Sommaire 


Modelling  of  U-tube  Tanks  for  ShipMoSD  Ship  Motion 
Predictions 

Kevin  McTaggart ;  DRDC  Atlantic  ECR  2011-300 ;  Recherche  et  developpement 
pour  la  defense  Canada  -  Atlantique ;  janvier  201 2. 

Introduction  :  Le  mouvement  de  roulis  d’un  navire  dans  les  vagues  pent  etre  impor¬ 
tant  en  raison  dn  faible  amortissement  et  de  la  proximite  de  la  frequence  naturelle  de 
roulis  du  navire  par  rapport  a  la  frequence  de  roulis  des  vagues  rencontrees.  Le  mou¬ 
vement  de  roulis  pent  etre  reduit  a  I’aide  de  divers  moyens,  y  compris  1’ installation 
de  quilles  de  roulis  et  de  dispositifs  de  stabilisation  actifs,  comme  des  gouvernails  on 
des  derives  de  stabilisation.  Les  citernes  a  tube  en  U  et  les  citernes  antiroulis  peuvent 
egalement  etre  utilisees  pour  reduire  le  mouvement  de  roulis ;  elles  ont  I’avantage 
d’etre  efficaces  meme  lorsque  le  navire  est  a  I’arret  et  sont  protegees  contre  les  dan¬ 
gers  associes  a  certaines  conditions  environnementales,  comme  la  glace. 

Resultats  principaux  :  Le  present  rapport  decrit  la  mise  en  oeuvre  d’un  modele  de 
citerne  a  tube  en  U  servant  a  enrichir  la  bibliotheque  de  donnees  sur  les  mouvements 
de  navires  du  logiciel  ShipMoSD.  Le  modele  de  citerne  a  tube  en  U  a  ete  mis  en 
oeuvre  pour  effectuer  des  calculs  dans  le  domaine  temporel  et  le  domaine  frequentiel. 
Des  exemples  de  calculs  pour  une  fregate  generique  demontrent  la  reduction  du  roulis 
lorsque  I’on  utilise  une  citerne  a  tube  en  U.  La  capacite  des  citernes  a  tube  en  U  a 
ete  integree  a  la  version  3.0  du  logiciel  ShipMoSD. 

Importance  des  resultats  :  La  bibliotheque  de  donnees  sur  les  mouvements  de 
navires  du  logiciel  ShipMoSD  pent  maintenant  etre  utilisee  pour  evaluer  I’incidence 
des  citernes  a  tube  en  U  sur  les  mouvements  de  navires.  Des  etudes  parametriques, 
y  compris  des  etudes  sur  I’incidence  des  dimensions  de  la  citerne  et  de  la  frequence 
naturelle  connexe  sur  le  mouvement  de  roulis  des  navires,  peuvent  etre  realisees.  Les 
previsions  dans  le  domaine  frequentiel  sont  efficaces  sur  le  plan  des  calculs  et  four- 
nissent  des  donnees  qui  sont  directement  utilisables.  Les  previsions  dans  le  domaine 
temporel  sont,  quant  a  elles,  utiles  pour  examiner  certains  phenomenes  non  lineaires, 
comme  c’est  le  cas  lorsque  le  mouvement  des  fluides  dans  les  citernes  a  tube  en  U  et 
I’efficacite  sont  limites  par  les  dimensions  de  la  citerne. 

Travaux  ulterieurs  prevus  :  La  bibliotheque  de  donnees  sur  les  mouvements  de 
navires  du  logiciel  ShipMoSD  sera  activement  mise  a  jour  et  appuyee. 
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1  Introduction 


Ship  motions  in  waves  affect  the  performance  of  both  humans  and  ship  systems.  Ship 
roll  motions  can  be  especially  large  because  the  natural  roll  frequency  of  a  ship  can 
be  similar  to  encountered  wave  frequencies  and  ship  roll  damping  can  be  quite  small. 
When  designing  a  ship,  several  steps  can  be  taken  to  reduce  roll  motions  in  waves. 
Bilge  keels  are  very  effective  and  are  widely  used.  Active  roll  stabilization  using 
rudders  or  dedicated  roll  stabilizer  hns  can  be  effective  if  a  ship  is  travelling  with 
sufficient  forward  speed  to  generate  required  lift  forces  on  rudders  or  hns.  U-tube 
tanks  and  hume  tanks  can  be  used  to  reduce  roll  motions,  and  are  effective  at  all  ship 
speeds.  Furthermore,  U-tube  tanks  and  hume  tanks  are  rugged  and  can  be  contained 
within  the  hull,  thus  making  them  protected  from  ice  and  other  environmental  hazards 
that  a  ship  might  encounter. 

In  anticipation  of  ongoing  design  work,  the  Canadian  Coast  Guard  tasked  DRDC 
Atlantic  with  adding  a  capability  for  modelling  U-tube  tanks  in  the  ShipMoSD  ship 
motion  library.  U-tube  tank  modelling  has  been  implemented  in  both  the  frequency 
and  time  domains  in  Version  3  of  ShipMoSD  [1,  2].  ShipMoSD  Version  3  also  includes 
modelling  of  sloshing  in  tanks  with  free  surfaces  [3].  The  U-tube  tank  implementa¬ 
tion  is  based  on  the  work  of  Lloyd  [4].  Section  2  gives  coordinate  systems  used  for 
evaluating  ship  motions  in  waves.  The  equations  of  motions  for  the  ship  and  U-tube 
tank  huid  are  given  in  Section  3,  with  the  numerical  implementation  described  in 
Section  4.  Section  5  gives  an  example  design  of  a  U-tube  tank  for  a  generic  frigate. 
Final  conclusions  are  given  in  Section  6. 

2  Coordinate  Systems  for  Ship  Motions  in 
Waves 


ShipMo3D  uses  2  main  coordinate  systems  for  evaluating  ship  motions  in  waves. 
Figure  1  shows  a  translating  earth  coordinate  system,  which  is  used  for  frequency 
domain  motions  and  also  for  evaluating  hydrodynamic  forces.  The  translating  earth 
coordinate  system  moves  with  the  mean  forward  speed  of  the  ship.  Figure  2  shows 
how  relative  sea  direction  is  dehned  in  the  translating  earth  coordinate  system.  For 
time  domain  computations,  which  can  involve  a  freely  maneuvering  ship,  motions 
are  given  in  the  earth-fixed  coordinate  system  of  Figure  3.  Within  ShipMo3D  time 
domain  computations,  forces  and  accelerations  at  each  time  step  are  evaluated  using  a 
translating  earth  coordinate  system  and  then  converted  to  the  earth-hxed  coordinate 
system. 
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Figure  1:  Translating  Earth  Coordinate  System  for  Frequency  Domain  Motions  and 
Solution  of  Hydrodynamic  Forces 


Figure  2:  Relative  Sea  Direction 
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Ship 


Figure  3:  Earth-Fixed  Coordinate  System  for  Time  Domain  Motions 

3  Equations  of  Motion  Inciuding  U-tube  Tank 


The  equations  of  ship  motion  including  a  U-tube  tank  can  be  evaluated  based  on  the 
work  of  Lloyd  [4] .  The  equations  presented  here  are  based  on  the  coordinate  systems 
presented  in  Figures  1  and  4;  thus,  some  terms  differ  from  those  in  Reference  4  due 
to  differences  in  coordinate  sign  conventions. 

For  a  ship  with  6  rigid-body  degrees  of  freedom  plus  a  degree  of  freedom  r  for  the 
U-tube  tank  (see  Figure  4),  the  equations  of  motion  in  the  frequency  domain  can  be 
written  as  follows: 


The  matrices  and  vectors  in  the  above  equation  are  partitioned  into  ship  motion 
and  U-tube  tank  fluid  angle  portions,  and  the  effective  mass  matrix  [M®]  includes 
contributions  from  both  ship  mass  and  added  mass.  The  following  term  is  introduced: 

Qtank  2  Ptank  '^r  ^  dj  tank  (2) 

where  ptank  is  the  density  of  fluid  in  the  tank,  Wr  is  the  reservoir  width,  w  is  the 
nominal  width  between  reservoirs,  and  Ltank  is  the  longitudinal  length  of  the  tank. 
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Figure  4:  Schematic  of  U-tube  Tank 
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Considering  first  only  the  displacement  of  the  tank  fluid,  the  effective  mass  and 
stiffness  terms  are: 


Mt 


a 


Qtank 


r 


Qtank  9 


(3) 

(4) 


where  hd  is  the  height  of  the  duct  and  is  the  nominal  height  of  the  fluid  in  each 
reservoir  relative  to  mid-height  of  the  duct.  The  undamped  natural  frequency  of  the 
tank  fluid  displacement  r  is  given  by: 


The  above  equation  can  be  re-expressed  using  the  tank  dimensions: 


(5) 


OOr 


ghd 


Wr  W  -\-  2  hr  hd 


(6) 


The  damping  associated  with  the  tank  fluid  displacement  r  is  highly  dependent 
on  the  local  fluid  flow  within  the  tank,  including  viscous  effects.  In  the  present 
formulation,  the  damping  Brr  is  expressed  in  terms  of  a  non-dimensional  damping 
coefficient  brr,  which  is  the  fraction  of  critical  damping: 


B 


TT 


2  b-fT  Qtank  9 

(jJr 


(7) 


The  damping  coefficient  b-rT  can  be  estimated  using  physical  model  tests  or  compu¬ 
tational  fluid  dynamics,  with  values  typically  being  of  the  order  of  0.10. 


The  ship  motions  influence  the  motion  of  the  tank  fluid,  giving  rise  to  the  following 
terms: 


Qtank 

(8) 

Qtank  (  ^r) 

(9) 

Qtank  ^tank 

(10) 

CrA 

Qtank  9 

(11) 

where  Zd  is  the  nominal  elevation  of  the  duct  relative  to  the  ship  centre  of  gravity  and 
Xtank  is  the  longitudinal  position  of  the  tank  relative  to  the  ship  centre  of  gravity. 

The  motion  of  the  tank  fluid  influences  the  ship  motions,  as  given  by  the  following: 


Qtank 

(12) 

Ml 

Qtank  (  “t“  ^r) 

(13) 

Ml 

Qtank  ^tank 

(14) 

CAr 

Qtank  9 

(15) 
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It  should  be  noted  that  the  ship  roll  stiffness  C44  does  not  include  the  effect  of 
the  U-tube  tank;  however,  the  effective  reduction  in  ship  metacentric  height  can  be 
evaluated  given  the  following  limit  for  low  frequency  motion: 

r  =  — ?74  for  cUe  — !■  0  (16) 

The  resulting  effective  change  in  static  metacentric  height  is 

A  ^  ^ 

where  A  is  the  ship  mass  displacement.  As  discussed  later,  the  above  effective  change 
in  metacentric  height  is  included  in  ShipMoSD  computations  with  a  U-tube  tank. 

The  fluid  motion  in  the  tank  will  be  limited  by  the  vertical  extent  of  the  tank,  which 
should  be  considered  during  both  design  and  time  domain  analysis.  The  maximum 
fluid  displacement  is: 

Tmai  =  arctan 

4  Implementation  within  ShipMoSD 


ShipMoSD  Version  3  has  been  completed,  including  user  manuals  [1,  2]  and  verifica¬ 
tion  and  validation  [5].  ShipMoSD  Version  3  introduces  the  following  capabilities: 

•  modelling  of  U-tube  tanks  as  described  in  this  report, 

•  modelling  of  sloshing  in  tanks  with  free  surfaces  [3], 

•  the  application  SMSDSeakeepSeawayFromRaos  for  fast  computation  of  ship 
motions  in  waves  using  previously  computed  response  amplitude  operators  [2]. 

Figure  5  shows  some  of  the  U-tube  tank  input  dimensions  for  the  application  SM3D- 
BuildShip  of  ShipMo3D  Version  3.  In  addition  to  modifications  to  SM3DBuildShip, 
the  ShipMo3D  applications  SM3DFreeMo,  SM3DSeakeepRegular,  SM3DSeakeepRan- 
dom,  and  SM3DSeakeepSeaway  have  been  modified  to  include  modelling  of  U-tube 
tanks. 

SM3DBuildShip  includes  an  input  record  for  correctionGM,  the  correction  to  meta¬ 
centric  height  GM.  This  term  can  be  used  to  account  for  effects  such  as  inter¬ 
nal  tanks  that  contain  liquid.  Note  that  the  input  for  a  U-tube  tank  in  Version 
3  of  SM3DBuildShip  accounts  for  the  dynamic  effects  of  the  U-tube  tank;  thus, 
correctionGM  should  not  include  effects  of  U-tube  tanks  for  which  full  input  has 
been  entered. 
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Figure  5:  ShipMoSD  Input  for  U-tube  Tank 
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Within  ShipMoSD  frequency  domain  computations,  the  solution  of  ship  motions  in¬ 
cluding  a  U-tube  tank  is  done  according  to  Equation  (1),  in  which  ship  and  U-tube 
tank  motions  are  considered  to  be  fully  coupled  and  are  solved  simultaneously.  When 
solving  for  ship  and  U-tube  tank  motions  in  the  time  domain,  minor  simplihcations 
are  used  such  that  the  ship  and  U-tube  tank  motions  are  solved  separately.  The 
following  is  assumed  when  solving  for  ship  motions  in  the  time  domain: 

f  {t  +  At)  f  (t)  (19) 

where  At  is  the  time  step  increment  for  solving  ship  motions.  The  above  assumption 
is  likely  acceptable  provided  that  At  is  small  or  the  inertia  forces  acting  on  the  ship 
due  to  the  U-tube  tank  are  small. 

The  motions  of  the  U-tube  tank  fluid  are  solved  in  the  time  domain  based  on  Equa¬ 
tion  (1): 

^  ^  (-  [Kr,]  {v}  -  [Brn]  {v}  -  Brrf  -  [Cr^]  {t]}  -  Crr  t)  (20) 

When  performing  time  domain  computations  of  the  fluid  displacement  angle  within 
a  U-tube  tank,  consideration  must  be  given  to  the  limiting  maximum  displacement 
angle  Tmax  imposed  by  tank  geometry.  If  the  limiting  maximum  displacement  angle 
will  be  exceeded  in  a  given  time  step  due  to  the  acceleration  given  by  Equation  (20), 
then  the  acceleration  is  modified  such  that  the  limit  displacement  angle  will  not  be 
noticeably  exceeded: 

^max  -  T{t)  -  f{t)At  ^  n 

=  (^tau - - for  T{t)  >  0  (21) 

-  for  r(t)  <  0  (22) 

where  etau  is  an  acceleration  attenuation  factor  for  maintaining  stability  during  time 
domain  simulations,  with  0.2  being  a  recommended  value. 

5  Computations  for  Generic  Frigate  with 
U-Tube  Tank 


Computations  have  been  performed  based  on  the  generic  frigate  of  the  ShipMoSD 
user  manuals  [1,  2].  These  sample  computations  demonstrate  the  design  of  a  U-tube 
tank,  and  also  provide  code  verification. 
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5.1  Design  of  U-tube  Tank 

Guidance  for  design  of  U-tube  tanks  is  provided  by  Lloyd  [4],  Gawad  et  al.  [6],  and 
Moaleji  and  Greig  [7].  For  optimal  dynamic  performance,  the  natural  frequency  of 
the  tank  fluid  (Equation  (5))  should  be  approximately  the  same  as  the  natural  roll 
frequency  of  the  ship.  To  limit  problems  with  loss  of  static  stability,  the  effective 
reduction  in  metacentric  height  (Equation  (17))  should  be  less  than  25  percent  of  the 
ship  metacentric  height.  Past  experience  indicates  that  a  U-tube  tank  will  have  a 
mass  of  1-5  percent  of  the  ship  mass.  Dynamic  performance  is  generally  best  if  the 
tank  is  at  a  similar  or  higher  elevation  than  the  ship  centre  of  gravity. 

Table  1  gives  the  properties  of  the  generic  frigate.  Using  these  properties,  the  design 
of  the  U-tube  tank  can  be  done  in  the  following  steps: 

1.  The  tank  is  located  at  midships.  The  bottom  of  the  tank  is  located  3  m  above 
the  baseline  (i.e.,  zu  =  3  m),  and  the  tank  has  a  total  height  ht  of  6  m.  The 
vertical  dimensions  of  the  tank  enable  large  fluid  motions  and  also  place  the 
tank  in  a  favourable  position  relative  to  the  centre  of  gravity.  The  total  width  of 
the  tank  Wd  +  ‘2wr  is  12  m,  taking  advantage  of  the  width  available  to  maximize 
effectiveness  for  available  volume.  The  reservoir  widths  Wr  are  assigned  values 
of  2  m,  leading  to  a  duct  width  tCrf  of  8  m  and  reservoir  spacing  w  of  10  m.  To 
maximize  the  allowable  fluid  displacement,  the  mean  fluid  level  h fluid  is  set  to 
{ht  -|-  hij)/2,  which  can  be  shown  mathematically  to  give  hr  =  0.5hf  (i.e.,  hr  = 
3  m). 

2.  The  effective  allowable  change  in  metacentric  height  is  specihed  to 

have  a  nominal  value  of  —0.2  m  (20  percent  of  a  projected  metacentric  height 
of  1.0  m  with  future  ship  growth).  The  parameter  Qtank  has  a  nominal  limit 
of  743000  kg  m  (Equation  (2)),  and  the  tank  length  Ltank  can  be  set  to  7  m 
(Equation  (2)). 

3.  The  tank  natural  roll  period  is  specihed  to  have  a  nominal  value  equal  to  the 
ship  natural  roll  frequency  of  0.704  rad/s.  Based  on  Equation  (6),  a  duct  height 
hd  of  0.6  m  can  be  used  to  give  a  tank  natural  frequency  of  0.706  rad/s. 

4.  The  mass  of  the  tank  huid  is  evaluated  to  be  126  tonnes,  3  percent  of  the  ship 
mass,  which  is  considered  to  be  acceptable. 

Table  2  gives  the  properties  of  the  example  tank,  and  Figure  6  shows  the  U-tube  tank 
and  the  hull  section  at  midships. 
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Table  1:  Generic  Frigate  Properties 


Length  between  perpendiculars,  Lpp 

120.0  m 

Draft  at  midships,  Tmid 

4.2  m 

Trim  by  stern,  tstem 

0.0  m 

Displacement,  A 

3713  tonnes 

Height  of  CG  above  baseline,  KG 

6.0  m 

Metacentric  height,  GM 

1.43  m 

Natural  roll  frequency,  UJ4 

0.704  rad/s 

Breadth  at  midships 

15.0  m 

Hull  depth  at  midships 

9.4  m 

Table  2:  Properties  of  Example  U-tube  Tank  for  Generic  Frigate 


Fluid  density  (fresh  water)  1000  kg/m^ 

Station  (AP  at  station  20)  10 

Length,  Ltank  7.0  m 

Duct  width,  Wd  8.0  m 

Reservoir  width,  Wr  2.0  m 

Height  of  bottom  above  baseline,  zu  3.0  m 

Total  height,  ht  6.0  m 

Duct  height,  hd  0.6  m 

Fluid  height,  hfiuid  3.3  m 

Fluid  mass,  m fluid  126  tonnes 

Tank  fluid  frequency,  ojtau  0.706  rad/s 

Maximum  fluid  angle,  T^ax  28.4° 

Effective  change  in  metacentric  height,  -0.189  m 
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Figure  6:  Example  U-tube  Tank  for  Generic  Frigate 


5.2  Verification  of  U-tube  Tank  Time  Domain 
Computations  with  Prescribed  Motions 

Numerical  tests  using  prescribed  motions  were  performed  to  verify  that  the  time 
domain  implementation  of  the  U-tube  tank  is  working  correctly.  All  tests  were  based 
on  the  example  U-tube  tank  with  properties  given  in  Table  2. 


A  decay  test  examined  fluid  displacement  angle  after  an  initial  fluid  displacement  of 
10  degrees.  The  U-tube  tank  itself  was  stationary  (i.e.,  not  subjected  to  ship  motion). 
Figure  7  shows  a  time  series  from  the  decay  test.  Many  standard  references,  including 
Reference  8,  give  equations  for  analyzing  dynamic  properties  based  on  time  series. 
The  damping  coefficient  based  on  the  time  series  was  computed  using: 

^  Iog(ry^°Vr^r) 

2  TV  n  UJ-r 


where  is  the  local  maximum  fluid  angle  after  i  cycles,  n  is  the  number  of  time 
series  cycles  for  analysis,  and  is  the  damped  natural  frequency  given  by: 


u 


damped 

tau 


Ulr  yi.O  - 


(24) 


The  natural  frequency  based  on  the  time  series  was  computed  using 


_  _ 2  TV  n _ 

'  ^  (‘“  -  C") 


(25) 


where  is  the  time  for  local  maximum  i. 
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Figure  7:  Time  Series  of  U-tube  Tank  Fluid  Angle  after  Initial  Displacement  Angle 
of  10  degrees 


Table  3  compares  the  actual  tank  natural  frequency  and  damping  coefficient  with 
values  obtained  from  analysis  of  the  decay  time  series.  The  values  obtained  from  the 
decay  time  series  are  in  excellent  agreement  with  the  actual  natural  frequency  and 
damping  coefficient  for  the  U-tube  tank. 

Table  3:  Comparison  of  U-tube  Tank  Natural  Frequency  and  Damping  with  Values 
Determined  by  Decay  Time  Series 


Natural  frequency  Damping  coefficient  b^-T- 
Actual  0.706  rad/s  0.100 

Decay  time  series  0.707  rad/s  0.100 


Two  additional  numerical  tests  considered  the  U-tube  tank  with  an  initial  fluid  dis¬ 
placement  angle  of  zero  and  the  tank  subjected  to  the  following  ship  roll  motion: 

riA{t)  =  0  for  f  <  0  (26) 

ri4{t)  =  fj4  [1.0  —  exp{t/tc)]  for  f  >  0  (27) 

where  is  the  roll  amplitude  and  tc  is  a  time  constant.  Figures  8  and  9  show  time 
series  for  ship  roll  and  fluid  displacement  angle  when  the  time  constant  tc  has  a  value 


12 


DRDC  Atlantic  ECR  201 1  -300 


of  10  s  and  the  roll  amplitude  has  values  of  15  degrees  and  40  degrees  respectively. 
Figure  8  shows  the  fluid  displacement  angle  approaches  a  value  of  —15  degrees 
as  expected.  In  Figure  9,  the  ultimate  fluid  displacement  angle  is  limited  by  the 
physical  limits  of  the  tank,  which  impose  a  limit  of  —28.4  degrees.  The  time  series 
for  the  large  roll  angle  indicates  that  the  time  domain  computations  can  manage  fluid 
displacement  angles  at  the  limits  of  the  tank. 

5.3  Verification  of  Ship  Motion  Computations 

Motion  computations  in  both  the  time  and  frequency  domains  were  performed  using 
the  generic  frigate  fitted  with  and  without  a  U-tube  tank.  Figures  10  and  11  show  roll 
RAOs  and  U-tube  tank  fluid  angle  RAOs  for  the  generic  frigate  travelling  at  10  knots 
in  beam  seas  with  wave  steepnesses  of  0.01  (1/100).  This  low  wave  steepness  was 
selected  so  that  the  U-tube  tank  fluid  motion  wouldn’t  be  limited  by  the  top  of  the 
tank,  thus  facilitating  comparison  between  frequency  and  time  domain  results.  The 
plotted  roll  motions  show  that  the  U-tube  tank  reduces  roll  amplitude  when  the  wave 
encounter  frequency  is  near  the  ship  natural  roll  frequency.  The  trends  in  Figure  10 
are  very  similar  to  those  presented  for  an  example  by  Lloyd  [4].  The  plotted  results 
show  excellent  agreement  between  computations  in  the  time  and  frequency  domains. 

Figures  12  and  13  show  RMS  roll  and  fluid  displacement  for  the  generic  frigate  trav¬ 
elling  at  10  knots  in  sea  state  5,  which  is  modelled  by  a  long-crested  seaway  and  a 
Bretschneider  spectrum  with  signihcant  wave  height  of  3.25  m  and  peak  wave  period 
of  9.7  s.  The  U-tube  tank  reduces  roll  motions  at  most  headings,  and  provides  the 
greatest  beneht  when  the  encountered  waves  have  a  period  near  the  natural  roll  fre¬ 
quency  of  the  ship.  It  should  be  noted  that  the  generic  frigate  has  substantial  bilge 
keels  to  provide  roll  damping.  For  a  ship  fitted  with  smaller  or  no  bilge  keels,  the 
relative  reduction  in  roll  motions  provided  by  the  U-tube  tank  would  be  significantly 
greater.  The  plotted  results  show  excellent  agreement  between  computations  in  the 
time  and  frequency  domains  for  the  ship  in  random  seas. 

5.4  Sensitivity  of  Motions  to  U-tube  Tank  Damping 
Coefficient 

When  modelling  a  U-tube  tank,  there  is  usually  uncertainty  when  selecting  a  value  for 
the  input  damping  coefficient  brT-  Figure  14  shows  frequency  domain  computations 
of  ship  roll  motions  using  tank  damping  coefficients  of  0.05,  0.10,  and  0.20.  The 
results  indicate  that  ship  roll  motions  are  not  very  sensitive  to  input  tank  damping 
coefficients  varying  within  a  realistic  range. 
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Tank  fluid  angle  r  (deg)  Roll  angle  r/4  (deg) 


Figure  8:  Time  Series  of  Roll  Motion  and  Fluid  Tank  Displacement  for  Ship  Roll 
?74  =  15°  [1.0  —  exp(— 1/10  s)] 
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Tank  fluid  angle  r  (deg)  Roll  angle  r  (deg) 


Figure  9:  Time  Series  of  Roll  Motion  and  Fluid  Tank  Displacement  for  Ship  Roll 
r/4  =  40°  [1.0  -  exp(-t/10  s)] 


DRDC  Atlantic  ECR  201 1  -300 


15 


Figure  10:  Roll  Amplitude  Versus  Wave  Frequency,  Generic  Frigate  at  10  knots  in 
Regular  Waves,  Beam  Seas  with  Steepness  of  0.01 


Figure  11:  U-tube  Tank  Fluid  Displacement  Amplitude  Versus  Wave  Frequency, 
Generic  Frigate  at  10  knots  in  Regular  Waves,  Beam  Seas  with  Steepness  of  0.01 
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Figure  12:  RMS  Roll  Versus  Relative  Sea  Direction,  Generic  Frigate  at  10  knots  in 
Sea  State  5 
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Figure  13:  RMS  U-tube  Tank  Fluid  Displacement  Angle  Versus  Sea  Direction, 
Generic  Frigate  at  10  knots  in  Sea  State  5 
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Figure  14:  RMS  Roll  Versus  Relative  Sea  Direction,  Generic  Frigate  at  10  knots  in 
Sea  State  5,  Different  U-tube  Tank  Damping  Coefficients 

5.5  Sensitivity  of  Ship  Motions  to  U-tube  Tank  Duct 
Height  and  Naturai  Frequency 

The  sensitivity  of  ship  motions  to  U-tnbe  tank  natural  frequency  has  been  examined 
by  varying  the  duct  height  for  the  tank  specihed  in  Table  2.  The  height  of  fluid 
h fluid  is  set  to  {hd  +  ht)  for  each  given  duct  height.  Table  4  gives  the  tank  natural 
frequencies  for  each  duct  height,  with  the  lowest  and  highest  natural  frequency  being 
within  20  percent  of  the  ship  roll  natural  frequency. 

Frequency  domain  computations  given  in  Table  15  indicate  that  the  roll  reduction  is 
generally  greatest  when  the  tank  natural  frequency  is  near  the  natural  roll  frequency 
of  the  ship. 


18 


DRDC  Atlantic  ECR  201 1  -300 


Table  4:  U-tube  Tank  Natural  Frequencies  for  Varying  Duct  Heights 


Duct  height  Natural  frequency 
bd  ^tau 

0.4  m  0.592  rad/s 

0.6  m  0.706  rad/s 

0.8  m  0.795  rad/s 


Figure  15:  RMS  Roll  Versus  Relative  Sea  Direction,  Generic  Frigate  at  10  knots  in 
Sea  State  5,  Different  U-tube  Tank  Duct  Heights 
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6  Conclusions 


Modelling  of  U-tube  tanks  has  been  implemented  into  the  ShipMoSD  ship  motion 
library.  Example  computations  for  a  generic  frigate  htted  with  a  U-tube  tank  exhibit 
expected  motion  behaviour  and  excellent  agreement  between  computations  in  the 
time  and  frequency  domains.  Uncertainty  typically  exists  regarding  selection  of  an 
input  damping  coefficient  for  a  U-tube  tank;  however,  example  computations  indicate 
that  ship  motions  aren’t  very  sensitive  to  variation  of  the  damping  coefficient  within 
a  realistic  range  of  values.  Example  computations  show  that  a  U-tube  tank  is  most 
effective  when  its  natural  frequency  is  near  the  ship  natural  roll  frequency. 
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Symbols  and  Abbreviations 


a  wave  amplitude 

[B]  damping  matrix 

brr  non-dimensional  damping  of  U-tube  tank 

6*^  non-dimensional  damping  of  U-tube  tank  from  time  series  analysis 

[C]  stiffness  matrix 

CG  centre  of  gravity 

g  gravitational  acceleration 

Hg  significant  wave  height 

hd  U-tube  tank  duct  height 

h fluid  height  of  fluid  in  U-tube  tank  relative  to  tank  bottom 

hr  nominal  fluid  height  in  U-tube  tank  reservoir 

ht  total  height  of  U-tube  tank 

k  wave  number 

Lpp  ship  length  between  perpendiculars 

Ltank  longitudinal  length  of  U-tube  tank 

[M®]  effective  mass  matrix 

m fluid  mass  of  U-tube  tank  fluid 

n  number  of  cycles  for  time  series  analysis 

Qtank  U-tube  tank  fluid  inertia  term 

RAO  response  amplitude  operator 

Tmid  ship  draft  at  midships 

tc  roll  motion  time  constant 

tstern  ship  trim  by  stern 

time  for  local  maximum  i 

w  nominal  distance  between  U-tube  tank  reservoirs 

Wr  U-tube  tank  reservoir  width 

X,  y,  z  translating  earth  coordinates 
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^tank 

Zbl 

Zd 

^-^tank 

At 


^tau 


Vj 

Va 

V 


Ptank 

T 

Tmax 

f 


X 

U4^ 


UJr 

damped 

^tau 


A 


earth-fixed  horizontal  plane  coordinates 
longitudinal  position  U-tube  tank  relative  ship  CG 
height  above  baseline 

nominal  elevation  of  U-tube  tank  duct  relative  to  ship  CG 
relative  sea  direction 

U-tube  tank  effective  influence  on  ship  metacentric  height 

time  step  for  solving  ship  motions 

U-tube  tank  acceleration  attenuation  factor 

motion  displacement  for  mode  j  in  translating-earth  coordinates 

roll  motion  amplitude 

wave  heading  (from) 

density  of  fluid  in  U-tube  tank 

U-tube  tank  fluid  displacement  angle 

U-tube  tank  maximum  possible  fluid  displacement  angle 

U-tube  tank  fluid  acceleration 

ship  heading  (to) 

ship  roll  natural  frequency 

natural  frequency  of  U-tube  tank  fluid  displacement 
damped  natural  frequency  of  U-tube  tank 
natural  frequency  of  U-tube  tank  from  time  series  analysis 
ship  mass  displacement 
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